Crossover from Non-Equilibrium to Equilibrium Behavior in the Time-Dependent 

Kondo Model 
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We investigate the equilibration of a Kondo model that is initially prepared in a non-equilibrium 
state towards its equilibrium behavior. Such initial non-equilibrium states can e.g. be realized in 
quantum dot experiments with time-dependent gate voltages. We evaluate the non-equilibrium 
spin-spin correlation function at the Toulouse point of the Kondo model exactly and analyze the 
crossover between non-equilibrium and equilibrium behavior as the non-equilibrium initial state 
evolves as a function of the waiting time for the first spin measurement. Using the flow equation 
method we extend these results to the experimentally relevant limit of small Kondo couplings. 



Introduction.- Equilibration in non-perturbative 
many-body problems is not well-understood with many 
fundamental questions still being unanswered. For 
example the crossover from a non equilibrium initial 
state to equilibrium behavior after a sufficiently long 
waiting time poses many interesting questions that are 
both experimentally relevant and theoretically of funda- 
mental importance. From the experimental side there is 
currrent interest in such questions related to transport 
experiments in quantum dots. Non-perturbative Kondo 
physics has been observed in quantum dots and given 
rise to a wealth of experimental and theoretical investi- 
gations. Quantum dot experiments allow the possibility 
to systematically study the effect of time-dependent 
parameters, like switching on the Kondo coupling at 
a specific time and measuring the time-dependent 
current. From the theoretical point of view this is 
related to studying the time-dependent buildup of the 
non-perturbative Kondo effect. 

The combination of strong-coupling behavior and 
time-dependent parameters makes such problems theo- 
retically very challenging. Various methods like time- 
dependent NCA and renormalized perturbation theory 
2] , the numerical renormalizationgroup [j| , bosonization 
and refermionization techniques Jj], etc. have allowed in- 
sights and e.g. identified the time scale t^ cx Ti/k^T^ 
related to the Kondo temperature 7k with the relevant 
time scale for the buildup of the Kondo resonance. Us- 
ing the form factor approach Lesage and Saleur could de- 
rive exact results for the spin expectation value P(t) = 
(S z (t)) for a product initial state pj- However, no exact 
results are available regarding the crossover from non- 
equilibrium to equilibrium behavior in this paradigm 
strong-coupling model of condensed matter theory. 

In this Letter we use bosonization and refermionization 
techniques to calculate the zero temperature spin-spin 
correlation function of the Kondo model at the Toulouse 
point |fj| exactly for two non-equilibrium preparations: 
I) The impurity spin is frozen for time t < 0. II) The 
impurity spin is decoupled from the Fermi sea for time 



t < 0. We find a crossover between non-equilibrium ex- 
ponential decay and equilibrium algebraic decay as one 
increases the waiting time for measuring the spin-spin 
correlation function at time t > 0. One concludes that 
zero temperature equilibration occurs exponentially fast 
with a time scale set by the inverse Kondo temperature, 
and a mixture of non-equilibrium and equilibrium be- 
havior for finite waiting times. Using the flow equation 
solution of the Kondo model pi |2| we then extend these 
results away from the Toulouse point to the experimen- 
tally relevant limit of small Kondo couplings in a system- 
atic approximation. 

Model. - The Kondo model describes the interaction of 
a spin-1/2 degree of freedom S with a Fermi sea 
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c 0a is the localized electron orbital at the im- 
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purity site. We allow for anisotropic couplings Ji — 
(J±,J±,J\\) and consider a linear dispersion relation 
(-k = vpk. Throughout this Letter we are interested in 
the universal behavior in the scaling limit Jj_/vp — > 0. 
The Kondo Hamiltonian can be mapped to the spin- 
boson model, which is the paradigm model of dissipative 
quantum mechanics 0- Our results can be interpreted 
in both these model settings; in the sequel we will focus 
on the Kondo model interpretation. 

We study two non-equilibrium preparations: I) The 
impurity spin is frozen for time t < by a large magnetic 
field term h(t)S z that is switched off at t = 0: h(t) ^> 7k 
for t < and h(t) = for t > 0. II) The impurity spin 
is decoupled from the bath degrees of freedom for time 
t < (like in situation I we assume (S z (t < 0)) = +1/2) 
and then the coupling is switched on at t = 0: Ji(t) = 
for t < and Ji(t) = Ji > time-independent for 
t > 0. This situation is of particular interest in fu- 
ture quantum dot experiments where the quantum dot 
is suddenly switched into the Kondo regime by applying 
a time-dependent voltage on a nearby gate [2| • 

The difference between these two initial states is that 
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in I) the electrons are in equilibrium with respect to the 
potential scattering induced by the frozen spin for t < 0. 
On the other hand in II) the initial state of the electrons 
is an unperturbed Fermi sea. We will later see that both 
initial states lead to the same spin dynamics. 

A suitable quantity for studying equilibration is 
the symmetrized zero temperature spin-spin correlation 

function C^ii(t w , t) d = \{{S z (t w ), S x (t w + In 
equilibrium this correlation function exhibits its well- 
known i _2 -algebraic long time decay and is, of course, 
independent from the initial (waiting) time t w : C eq (t) = 
C e q(t w , t) oc t 2 . 

Exact results for the non-equilibrium spin dynam- 
ics have so far only been obtained for the spin expec- 
tation value P(t) = f (SJt)). At the Toulouse point 
J^ p /2nv F = 1 - 1/V2 ja] P(t) can be evaluated ex- 
actly and one finds a purely exponential decay P(t) = 
h exp(— 2t / ■Kwt^) . Here and in the sequel, we define the 
Kondo time scale as £k = 1 /2k with the Kondo temper- 
ature defined via the impurity contribution to the Som- 
merfeld coefficient 7i mp = wit 2 /3Tk, where w = 0.4128 
is the Wilson number. Using the form factor approach 
Lesage and Saleur could derive exact results for P(t) even 
away from the Toulouse point For < J\\ < Jl F 
they find that the spin expectation value again decays 
exponentially for large times t/tx 3> 1 with the same ex- 
ponential dependence P(t) oc exp(— 2t/iTwt-^), however, 
its behavior at finite times is more complicated. 

Since P(t) = 2C(t w — 0, t) these results raise the ques- 
tion of how the crossover between non-equilibrium expo- 
nential decay and equilibrium algebraic decay occurs as 
a function of the waiting time t w . 

Toulouse point. - We have addressed this issue by eval- 
uating Ci.u(t w ,t) exactly at the Toulouse point of the 
model. One finds that the result is the same in both 
situation I) and II): 

Ci t ii(t w ,t) — C eq (t) 

-2e- t "'/* B s(t)(s(^)e- t/tB - s{t w + t)) 
_ e -2W*B ( s (f w ) e -*/tB _ s{tw + t) y (2) 

for t,t w > 0. Here s(r) = (£b/tt) J °° dw sin(W)/(l + 

w 2 ig) with the abbreviation £b = f 7ru;i K . In terms of 
s(t) the equilibrium correlation function reads 

C^{t) = \e-^-s\t) . (3) 

Notice that sit) ~ t&/wt for t ^> t-Q leading to the alge- 
braic long-time decay. Therefore the Fourier transform 
of the equilibrium correlation function is proportional to 
\oj\ for small frequencies uj/Tk <C 1: C eq (w) oc |w|/T|-. 

For zero waiting time t w = the non-equilibrium cor- 
relation function Eq. shows the well-known purely 
exponential decay Ci t n(t w — 0,t) — je~ 2t ^ tB , while for 
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FIG. 1: Universal curves for the spin-spin correlation func- 
tion C(t w ,u>) at the Toulouse point for various waiting times 
(t m — 0, tyi/A, tx/2, tn, 2tK, 00 from top to bottom). The inset 
depicts the same curves on a linear scale, which shows more 
clearly the onset of the non-analyticity (C(t m ,u>) — C(t m , 0) oc 
|w|) for waiting times t w > leading to the algebraic long- 
time decay. 

any fixed t w > the algebraic long time behavior domi- 
nates with an amplitude that is suppressed depending on 
the waiting time: 

d, u (t w , t) ~ ~s 2 (t) (l - e -Wta) 2 + o ((tn/t) 3 ) ( 4 ) 

for t 3> xxxax(t w , is)- In particular, one can read off 
from 121 that the difference between the non-equilibrium 
and equilibrium correlation function decays exponen- 
tially fast as a function of the waiting time 

|Ci,ii(W) - Ceq(*)| oce-*»/* B (5) 

for t w 3> ts • These results can be understood by noticing 
that the initial state corresponds to an excited state of 
the model, therefore yielding a different spin dynamics 
from equilibrium. After a time scale of order £k corre- 
sponding to the low-energy scale of the model, the ini- 
tial non-equilibrium state leads to the behavior of the 
equilibrium ground state with deviations that decay ex- 
ponentially fast as one increases the waiting time |ld| . 
This crossover behavior of the non-equilibrium correla- 
tion function at the Toulouse point is shown in Fig. Q] for 
the one-sided Fourier transform with respect to the time 
difference t 

f°° dt 

C(t w ,u;)= ^ e ^ C(t w ,\t\) . (6) 

Method.- We perform the standard procedure of 
bosonizing the Kondo Hamiltonian in terms of spin- 
density excitations b k , b\ , and then eliminating its 
Jii-coupling using a polaron transformation U — 

exp[i<5 S z J2k>o( 1 /^) ( b k ~ b l)} ( for details see e.g. 0). 
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At the Toulouse point the transformed Hamiltonian H = 

, r— l 

W HU can be refermionized 9] 

H = J2^ k + Y, v ^{ xS 'l d+d ^ k ) ■ ( ? ) 

k k 

Here the spinless fermions ^(x) correspond to soli- 
ton excitations built from the bosonic spin-density 
waves with the refcrmionization identity ^(x) oc 

exp [iEkyoi 1 /^) (K e ~ lkx ~ ^e 4 **)]. Eq. can be 
interpreted as a resonant level model with the hybridiza- 
tion function A(e) = f ^J2k^k^( e ~ e k) = Tk/kw and 
a fermionic impurity orbital d with the identity S z = 
d)d — 1/2. Model I is therefore represented by a reso- 
nant level model with a time-dependent impurity orbital 
energy e d (t) = h(t) 

Hj(t) = 5> fc *l* k + <-,,(*) (dtd- 1/2) 

k 

+ Y / Vk(^ k d + d^ k ) (8) 

k 

with td{t < 0) = oo, ed(t > 0) = 0. In model II the spin 
is not coupled to the Fermi sea for time t < 0: this leads 
to a time-dependent hybridization function and a time- 
dependent potential scattering term (due to the polaron 
transformation U) 

Hu(t) = £ e fc *l* fc + £ Vk (t) (*t d + rft* fc ) 
fe fc 

+ (5 2 (i = 0))^ 5fcfc ,(i)(*I* fe '-l/2) (9) 
fc.fc' 

with < 0) = (1 - t/2)/2wvf, Qkk'(t > 0) = and 

A(e; i < 0) = 0, A(e; t > 0) = T K /W 

In order to evaluate the non-equilibrium spin dynamics 
we use the quadratic form of ffi,n to solve the Heisenberg 
equations of motion for the operators cft(t),d(t) exactly. 
After some straightforward algebra one can write the op- 
erator S z (t) for t > as a quadratic expression in terms 
of the operators ^\(t = 0), * fe (f = 0), S{t = 0) and 
d(t = 0) at time t = |ll| . The non-equilibrium corre- 
lation function is then given by 

Cui(t w ,t) = ^(G 1 j 1 \{S z (t w ),S z (t w + t)}\G 1J1 ) . (10) 

where we insert the solution of the Heisenberg equation 
of motion for S z (t). The initial non-equilibrium states 
|Giji) remain time-independent in the Heisenberg pic- 
ture and are simply the ground states of the Hamiltoni- 
ans Hiji(t) for time t < 0. Inserting these expressions 
m (fTTTfl yields 101 after some tedious but straightforward 
algebra . For completeness we also give the result for 
the imaginary part of the non-equilibrium Greens func- 
tion 

K ul (t w , t) d = 1 ( [S z (t w ), S z (t w + t)} ) U i (11) 
= -ze-*/* B (s(t) - s(t w + t)e-*™/* B + s(t w )e-^+W tB ) 



which again approaches the equilibrium result K eq (t) = 
—i e~ t t tB s{t) exponentially fast as a function of t w /t^. 

Kondo limit.- The Toulouse point exhibits many uni- 
versal features of the strong-coupling phase of the Kondo 
model like local Fermi liquid properties, however, other 
universal properties like the Wilson ratio depend explic- 
itly on the coupling Ju . This raises the question which of 
the above non-equilibrium to equilibrium crossover prop- 
erties are generic in the strong-coupling phase. We inves- 
tigate this question by using the flow equation method 
[l2j that allows us to extend our analysis away from the 
Toulouse point in a controlled expansion. In this Let- 
ter we focus on the experimentally most relevant limit 
of small Kondo couplings (notice that the flow equation 
approach is not restricted to this limit and can be used 
for general J\\ 0). 

The flow equation method diagonalizes a many- 
particle Hamiltonian through a sequence of infinitesimal 
unitary transformations in a systematic approximation 
[l2^ . This approach was carried through for the Kondo 
model in Ref. 0. Since the Hamiltonian is transformed 
into its diagonal basis, we can follow the same steps as 
in the Toulouse point analysis: i) The Heisenberg equa- 
tions of motion for the unitarily transformed observables 
can be solved easily with respect to the diagonal Hamil- 
tonian. ii) One then re-expresses the time-evolved oper- 
ators through the operators in the initial (non-diagonal) 
basis for time t = 0. iii) The correlation functions (|1C)|) 
for general t w are evaluated. An operator product ex- 
pansion to leading order is employed like in Ref. to 
close the resulting systems of equations. 

The above procedure can be used quite generally to ap- 
ply the flow equation method to time-dependent Hamil- 
tonians of the above kind. For the Kondo model specifi- 
cally, one can simplify the calculation by using the results 
from Ref. 0: it was shown that a resonant level model 
(|Tj) with a universal non-trivial hybridization function 
Aeff(e) 7^ const, can be used as an effective model for the 
spin dynamics on all time scales; the only free parame- 
ter is the low-energy scale Tk- Similarly, the effective 
Hamiltonian © with A(e;t > 0) = A cff (e) from Ref. 
yields the .S^-spin dynamics for both non-equilibrium 
situations I and II. A careful analysis 0| shows that the 
only effect not captured by the resonant level model is 
the polaron-like transformation that is contained in the 
complete flow equation approach. This polaron tranfor- 
mation leads to an initial potential scattering term like 
in with gkk'(t < 0) = (A(-B e ff) - V2)/2itvf, where 
B c g = l/( e fc + e l') an d A(£Qis the flowing scaling di- 
mension according to Ref. [7j = 0) = \/2 and 
X(B — ► 00) = 1). However, this initial potential scat- 
tering term has a neglible effect (relative error < 5%) 
unless one is interested in short waiting times t w < £k/4: 
for larger waiting times the effective Hamiltonian from 
Ref. |8| yields a very accurate description without the 
need for the full flow equation analysis. 
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FIG. 2: Universal curves for the spin-spin correlation func- 
tion C(t w ,uo) in the limit of small Kondo couplings (Kondo 
limit) for various waiting times (t w — 0, tjc/4, ix/2, ix, 2£k, oo 
from top to bottom like in Fig. . The inset shows the spin 
expectation value P(t): the dashed line is an asymptotic fit 

^asym 

(t) = 0.11exp(-1.51t/t K ). 

From the quadratic Hamiltonian @ with the above 
effective coupling constants (including the initial poten- 
tial scattering term) one can easily evaluate the non- 
equilibrium correlation functions; results are depicted 
in Fig. |21 The key observations from the Toulouse 
point analysis hold in the Kondo limit as well, only the 
crossover behavior is more complicated: 1) The system 
approaches equilibrium behavior exponentially fast as a 
function of t w /tK (compare Eq. (J5J). Notice that the ini- 
tial approach for small t w in Fig. [21 is faster than at the 
Toulouse point (Fig. QJ. 2) An algebraic long-time de- 
cay oc t~ 2 dominates for all nonzero waiting times t w > 
(compare Eq. Q). 

For zero waiting time t w = the inset in Fig. [21 
shows the decay of the spin expectation value P(t) in 
the Kondo limit, which has not been previously calcu- 
lated explicitly on all time scales. For large t/tK the 
behavior crosses over into an exponential decay which 
agrees very well with the exact asymptotic result from 
Ref. 0: P C xact(t) oc exp(— 2t/-Kwt-K) oc exp(— l.bAt/tyi). 
On shorter time scales the decay is faster, which is due to 
unrenormalized coupling constants at large energies that 
dominate the short-time behavior. 

Conclusions.- Summing up, we have investigated the 
crossover to equilibrium behavior for a Kondo model that 
is prepared in an initial non-equilibrium state. We cal- 
culated the non-equilibrium spin-spin correlation func- 
tion on all time scales and could show that it evolves 
exponentially fast towards its equilibrium form for large 
waiting time of the first spin measurement t w 3> £k, see 
Eq. (JjjJ. Our results also established that the flow equa- 
tion method is a very suitable approach for studying such 
non-equilibrium problems: it agrees with very good ac- 
curacy with exact results for both t w = and 2™ = oo, 



and can describe the crossover regime as well. 

Finally, it is worthwile to recall the fundamental quan- 
tum mechanical observation that the overlap between the 
time-evolved non-equilibrium state and the true ground 
state of the Kondo model is always time-independent. 
Therefore it is not strictly accurate to conclude from 
our results that an initial non-equilibrium state "de- 
cays" into the ground state: rather, quantum observables 
which exhibit equilibration behavior are probes for which 
the time-evolved initial non-equilibrium state eventually 
"looks like" the ground state. Since the notion of equili- 
bration into the ground state plays a fundamental role in 
quantum physics, it would be very interesting to study 
other systems with quantum dissipation to see which of 
the crossover and equilibration properties derived in this 
Letter are generic. 
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